From 3D to 2D Hydrodynamics in Interacting Micro-rods 
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Moving micron scale objects are strongly coupled to each other by hydrodynamic interactions. 
The strength of this coupling decays as the inverse particle separation when the two objects are 
sufficiently far apart. It has been recently demonstrated that the reduced dimensionality of thin 
fluid layer gives rise to longer ranged, logarithmic coupling. Using holographic tweezers we show 
that microrods display both behaviors interacting like point particle in 3D at large distance and 
like point particles in 2D for distances shorter then their length. We derive a simple analytical 
expression that fits remarkably well our data and further validate it with finite element analysis. 



Flagella, microtubules, nanotubes and nanowires are 
a few examples of the rich variety of fundamental roles 
that slender bodies have in physics, chemistry and biol- 
ogy. They all share the same geometrical feature of be- 
ing "slender" , that is having a linear dimension which is 
much larger than the other two. As a consequence of that 
they're all expected to have a similar physical behavior 
in those situations where only shape matters. That is, 
for example, the case for micro- hydrodynamics. A thin 
nanotube or a thicker microrod will both experience a 
similar drag force, which will depend mainly on the body 
length while thickness only appears in logarithmic terms 
[T]. While there's a huge amount of work on single slen- 
der body dynamics, especially in the context of bacterial 
motility [2], investigation of hydrodynamic interactions 
between slender bodies is almost only limited to the phe- 
nomenology of syncronization [3l [4] . One reason could 
be that, as opposed to spheres, which have been stud- 
ied extensively [SHZ], hydrodynamic couplings between 
anisotropic bodies are a complex function of both rel- 
ative distance and orientation. Optical tweezers can be 
used to trap and move one dimensional objects [SHTO], In 
particular, holographic optical trapping has been shown 
to be an ideal tool for full 3D micromanipulation of mi- 
crorods and nanotubes pTUTS] , Such capabilities offer 
a unique opportunity for trapping and orienting slender 
bodies in well defined relative configuration, and directly 
probe their coupled Brownian dynamics. 

In this Letter we provide a direct measurement of hy- 
drodynamic coupling between a pair of parallel aligned 
silica microrods optically trapped in blinking holographic 
tweezers. We found a crossover from 3D behavior at large 
distances, to a 2D logarithmic behavior, when the dis- 
tance falls below the rods length. Experimental data are 
in excellent agreement with finite element analysis and 
can be very well reproduced by a simplified theoretical 
approach. 

In a practically zero Reynolds number regime, as it 
is in the mesoscopic world, a moving objects generates 
a perturbation in the surrounding fiuid that decays as 
the inverse distance (TJ [16]. A nearby particle will be 



transported by this fiow and experience a higher mobil- 
ity when trying to move along the same direction of the 
other particle. For the same reasons relative motions 
are characterized by lower mobilities. For well separated 
objects, collective and relative mobilities can be easily 
derived assuming that each object is rigidly transported 
by the approximately uniform flow field produced by the 
other one. Within this approximation when forces Fi 
and F2 are applied to two bodies they will move with 
speeds: 



Ui = Ml -Fi +U12 
U2 = M2 • F2 + U21 



(1) 



where is the ith body mobility tensor while Uij rep- 
resents the fiow that object j produces at the location 
of i. Due to linearity in Stokes fiow, flow fields have the 
form: 



(2) 



the fiow propagator Gij will depend on distance and 
on the relative orientation of particle j. Introducing 
the velocity and force vectors U = (Ui,U2,...) and 
F = (Fi,F2,...), the full many-body problem can be 
then stated in the form of a compact mobility matrix 
formulation: 



U = ]M-F 



where M is the mobility matrix: 
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with hydrodynamic couplings appearing as oflF-diagonal 
terms. At large enough distances, whatever is the shape 
of j, the propagator will tend to to the Oseen tensor that 
only depends on the position vector Vij = — r : 
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The Oseen tensor represents the first term in the mul- 
tipole expansion of the Stokes flow produced by a given 
force spatial distribution, also known as the Stokeslet. 
The Stokeslet propagates the perturbation of a zero- 
dimensional point force and therefore does not depend 
on particle shape and orientation. For spherical beads 
it provides a remarkably good description of couplings 
down to interparticle distances of about 2.5 radii [SHZ]- 
For anisotropic bodies, higher order terms will have to 
be included and their form will depend on particle shape 
and orientation. A thin rod may be thought as the sim- 
plest anisotropic body consisting of a one dimensional 
line force distribution. Most of the shape dependence 
occurs via a single parameter, that is the body length L, 
and rod orientation matters. Such a dependence on ori- 
entations, coupled to the long range character of Stokes 
propagators, makes hydrodynamic coupling between thin 
rods a very challenging problem both from the experi- 
mental and theoretical point of view. As a first start- 
ing approach we could limit ourselves to a single pair of 
interacting rods at a fixed distance. Slender body the- 



ory provides a practical theoretical framework for uni- 
dimensional objects where the point force description of 
the sphere is substituted by a line distribution of point 
forces. Slender body theory has been successfully applied 
in the context of flagellar and ciliar propulsion in bacte- 
ria. We have seen that coupling terms in the mobility ma- 
trix can be thought, as a first approximation, as the flow 
propagators between interacting objects. Although this 
is straightforward in the case of spheres, where sphere's 
center is a good representative point of sphere position 
and center of forces, it is quite ambiguous when bodies of 
extended size are considered. It can be shown [16] that 
when a thin prolate spheroid is immersed in an exter- 
nal inhomogeneous flow, it will move with a speed given 
by the average of the original flow over the rod's length. 
This finding suggests that we can approximate the hydro- 
dynamic coupling between the rods as the average flow 
produced by one rod over the length of the other one. 
For parallel aligned rods like in Fig. [l] this amounts to 
the double Stokeslet integration: 
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FIG. 1: a) Geometrical parameters of the investigated micro- 
rods configurations, b) Two silica microrods, each one held 
in a double holographic trap. 

When the two rods are much far apart than their 
length, the obtained propagator reduces obviously to the 
Oseen tensor component: 
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than rods length, the propagator has the logarithmic ex- 
pression: 
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More interestingly, for separations that are much smaller 



Such a form reminds very closely of what is observed 
for hydrodynamic coupling in two dimensional thin liq- 
uid layers [17^. In other words, microrods interact like 
point particles in 3D when they are very far away, and 
then like point particles in 2D at short distances. The 
crossover distance is determined by the most important 
length in the problem, that is the rod's length. Such a 
finding is not surprising when one realizes that at short 
distances the problem reduces to the idealized two di- 
mensional case of infinitely long cylinders, where the fi- 
nite length of the two rods is only responsible for edge 
effects. Our simple formula seems to provide a good in- 
terpolation between the two limiting cases but too crude 
approximations might be involved in its derivation and 
we decided to test it against finite element analysis. The 
model consists of two cylinders of length L and radius a 
at distance d. Finite size effects are reduced by having 
the two cylinders inside a fluid sphere of radius lOOL. 
No-slip boundary conditions are set on the cylinder sur- 
faces. We then calculate Stokes drag on the cylinders in 
the case of rigid and relative motion at different d values. 
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Hydro dynamic coupling is then obtained as half the dif- 
ference between rigid and relative mobilities. The result- 
ing values for coupling G^^ are reported as solid circles 
in Fig. [2] showing a remarkable good agreement with the 
analytical expression in ([6| also plotted as a solid line. 

The chance of observing two rods, occasionally aligned 
in parallel at some given distance, clearly makes the prob- 
lem practically impossible to approach experimentally 
by simple video microscopy observations. However holo- 
graphic tweezers provide an ideal tool to directly verify 
our results by trapping and orienting two slender objects 
in well defined and reproducible relative configurations. 
We used holographic tweezers to manipulate silica mi- 
crorods of about 10 jiia in length and 150 nm in radius 
[15]. Multiple optical traps are obtained by shaping an 
expanded CW laser beam (532 nm) with a spatial light 
modulator (Holoeye LC-R 2500). The shaped beam is 
then focused by a microscope objective of high numer- 
ical aperture (Nikon Plan Apo lOOx, NA 1.4) onto the 
target array of trapping spots. Optimal intensity distri- 
bution are generated by GSW algorithm running on an 
NVIDIA GPU (GTX 480) HHH]. Each rod is held in 
two traps independently and dynamically reconfigurable 
(Fig. ^p). The two rods can be aligned and moved 
at different distances. We limit ourselves to the com- 
puted case of center of mass dynamics of parallel aligned 
rods. Using a chopper on the trapping beam, we peri- 
odically release the rods every 1/15 s and subsequently 
start recording bright field images at 300 fps for 1/30 
s. The horizontal x coordinate of each rod is tracked by 
fitting the row average intensity profile of cropped rods 
images. A straightforward way to obtain hydrodynamic 
couplings is that of tracking rods free Brownian motion. 
We have estimated and experimentally verified that our 
rods have a rotational diffusion coefficient of about 0.02 
rad^/s. That means that the mean squared angular dis- 
placement during acquisition time is about tt/IOO. It is 
therefore safe to neglect rotational-translational coupling 
[2m I2T] and project on the x axis the mobility problem 
with stochastic forces: 
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where mi and m2 are the transverse mobilities of the 
two rods while r]i and 772 are stochastic forces with zero 
average and correlation: 



(10) 



This is a coupled Langevin equation leading to a dis- 
placements covariance matrix [22 : 
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(11) 




FIG. 2: Hydrodynamic coupling G^^ between two parallel 
aligned thin rods as a function of their relative distance. Solid 
line represents the simple theory in eq. (|6|. Black circles are 
numerical results from finite element analysis. Experimental 
data are reported as open squares. Dashed lines represent 
the long ^ and short Q distance limiting behavior reducing 
to respectively 3D and 2D hydrodynamic couplings between 
point particles. 



Diagonal terms represent the mean squared displace- 
ments of the two rods. The corresponding diffusion co- 
efficients provide a direct measurement of single particle 
mobilities. On the other hand hydrodynamic couplings 
are easily extracted from the diffusion coefficients of the 
off diagonal terms. In Fig. [3] we report the two sin- 
gle particle mean squared displacements toghether with 
the crossed term for interparticle distances of 3.5 and 20 
microns. Diffusion coefficients can be easily extracted 
by fitting with a straight line. Single particle diffusion 
coefficients don't show a systematic dependence on in- 
terparticle distance and their deviations are mainly at- 
tributable to small differences in rod lengths. Microrods 
in our sample have an average length of 10.5 /im with a 
standard deviation of 2.5%. Tracking single rod Brown- 
ian motion, we extract an average transverse mobility of 
42 fim/s pN. By theoretical predictions, this corresponds 
to an average rod thickness of about 200 nm, which is well 
compatible with the nominal pore size (300 nm) used for 
rod's growth [15]. Turning now to hydrodynamic cou- 
plings, the diffusion coefficient of crossed terms directly 
provides a measure of G^^{d) through (11). Indeed we 



found an excellent agreement with theoretical and numer- 
ical predictions as shown in Fig. |2] It is worth to note 
that coupling values are expected to be much less affected 
by the actual rod thickness than the single rod mobilities, 
as can be understood noticing that the thickness never 
entered in the arguments leading to (|6|. In particular, 
the results in Fig. [2] are expected to remain valid even 
in the case of vanishing rods thickness, as could be the 
case for single walled carbon nanotubes, microtubules or 
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short straight DNA segments. 
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FIG. 3: Mean squared displacements of freely diffusing rods 
starting off a parallel aligned configuration generated with 
holographic optical tweezers. Top and bottom panel refer re- 
spectively to interparticle distances of 3.5 and 20 /xm. Solid 
circles and squares represent the diagonal terms of the gen- 
eralized MSD while solid triangles refer to the off diagonal 
terms. Linear fits are also reported as solid lines. 

In conclusion, we have directly measured hydrody- 
namic interactions between freely diffusing micro-rods. 
We demonstrated by direct experimental observation and 
numerical finite element analysis that two parallel aligned 
microrods interact like point particles in 3D when they're 
much farther then their lentgh, and like point particles in 
2D for short distances. We also derived a simple analyti- 
cal expression that reproduces very well both experimen- 
tal and numerical data. Although the validity of our re- 
sults is expected to hold even for much thinner, and inter- 



esting objects, like nanotubes or microtubules, an exper- 
imental investigation of hydrodynamic coupling in those 
systems remains an open and challenging problem. This 
work was partially supported by IIT-SEED BACTMO- 
BIL project and MIUR-FIRB project RBFR08WDBE. 
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